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Abstract 

We introduce a fractional Fokker-Planck equation with a temporal power-law dependence on 
the drift force fields. For this case, the moments of the tracer from the force-force correlation 
in terms of the time-dependent drift force fields are discussed analytically. The long-time 
asymptotic behavior of the second moment is determined by the scaling exponent £ imposed 
by the drift force fields. In the special case of the space scaling value v = 1 and the time 
scaling value r = 1, our result can be classified according to the temporal scaling of the mean 



second moment of the tracer for large t: < x 2 {t) > oc t with £ = 7 for normal diffusion, and 



< x 2 {t) > oc t v with 77 > 1 and £ > 7 for superdiffusion. 
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Recently, a lot of interest has been concentrated on nonstationary random processes in 
the behavior of disordered systems [1 — 3]. The anomalous feature of diffusive stochastic 
processes is known to be qualitatively different from the standard behavior characteristic of 
the regular system. Especially the anomalous diffusive behavior in the disordered system 
has mainly been argued from the random walk, a useful theory of classical stochastic pro- 
cesses in the past; in fact, the normal diffusion of tracers is well known to be a result of 
random walkers characterized by the variance < x 2 {t) > oc t. The random walk theory with 
the nearest-neighboring transition probability has been numerically discussed by using the 
Green function method. The anomalous transport process of the random walk theory has 
been widely extended to a continuous-time random walk theory, which is specified by the 
distributions of both the pausing time and the transition probability, which are dependent on 
the length between steps. It provides a dynamic framework for a precise connection between 
fractional diffusion equations and fractal walks [4, 5] while the transport phenomena for the 
motions of tracers have largely been extended to reaction kinetics[6, 7] and strange kinetics 
[8]. 

It is well known that the second moment for the super diffusive proess is proportional to 
t a with a > 1 for superdiffusion and a < 1 for subdiffusion and to log^t for strong diffusion. 
In the past, the best example of the super diffusive random proess was a classical paper by 
Richardson [9] in which the second moment for a relative separation I at time t, < I 2 > oc 
t 3 , was shown to depend on the motion of two particle in a turbulent fluid. Matheron and 
de Marsily [10] showed the anomalous behavior of the longitudinal dispersion, < x 2 (t) > 
oc t 3 / 2 , in a stratified random velocity field. Among other recent examples, we also can 
mention anomalous kinetics in chaotic dynamics due to flights and trapping [8], layered 
velocity fields [11], anomalous diffusion of amphiphilic molecules [12], anomalous diffusion in 
a two-dimensional motion of rotating flow [13], etc. Examples for subdiffusive stochastic 
processes are charge transport in amorphous semiconductors [14, 15], NMR diffusometry on 
percolation structures [16], and the dynamics of a bead in polymers [17]. Furthermore, the 
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stochastic process for anomalous diffusion on a fractal structure and Levy freights has been 
the subject of several reports in the literature [18—20]. The typical models using Levy heights 
include bulk-mediated surface diffusion [21], transport inmicelle systems or heterogeneous 
rocks [22], single molecule spectroscopy [23], and special reaction dynamics [24]. Among the 
many frameworks for characterizing anomalous diffusion are fractional Brownian motion [25], 
fractional diffusion equations [26, 27], and generalized Langevin and Fokker-Planck equations 
[28-31]. 

The motivation for this paper was to apply the analytical method of a field-theoretical 
space-time renormalization group to investigate the anomalous behavior of the fractional 
diffusion equation. The moments for the probability density will be found so that the prop- 
erties lead to anomalous asymptotics of the diffusion process. The purpose of this paper 
is to present an anomalous transport process for a test particle described by a fractional 
Fokker-Planck equation in the presence of a temporal power-law dependence for the force- 
force correlation of the drift force fields. The resulting formulae for the diffusive behavior of 
the tracer in the Fourier-Laplace region are discussed analytically, and the anomalously long- 
time behavior of the moments can be determined by using the scaling exponent £ imposed 
by drift force fields. 

First, the fractional diffusion equation for the probability density P(r,t) on the d- 
dimensional fractal structure can be expressed in terms of 

d 13 1 

W P{r,t) = ^Vai^Or^VaPfcf), (1) 

where V a = d ®^ a and D(r) is the anomalous diffusivity and depends only on the position. 
The normalization condition for the probability density is given by 

poo 

/ Dr d - 1 P(r,t)dr = l. (2) 
Jo 

If a — (3 — 1, the form of Eq. (1) transforms to a spherically symmetric diffusion equation 
in (i-dimensional Euclidean space, showing that this equation has analytical solutions for 
a diffusive behavior characteristic of fractal systems. Let us now introduce the anomalous 
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diffusivity given by 



D(r) = Dr\ 



(3) 



where D is the diffusion constant and 9 is the scaling exponent and can be determined from 
the fractal structure. The feature of this diffusivity describes the statistical property of 
diffusion on fractals, and even though it extends to a number of 9 and a we here consider 
the special case of 9 = a and < a < 1 in order to find the moments. After multiplying Eq. 
(1) by r a and r 2a , respectively, and integrating, we obtain the asymptotic result 



This represents the behavior for a super diffusive particle for the scaling exponent relation 
2(3 > a and a subdiffusive particle for < 2(3 < a. From Eq. (4), the value of the scaling 
exponent, 2(3 = a is equal to one for normal diffusion. 

From now, we introduce the equation of motion in the presence of a spatiotemporal 
power-law dependence of both the drift force and an external harmonic potential. Since 
we take into account an overdamped particle of mass m moving in a viscous medium. The 
equation of motion can be expressed in terms of 



where f(v, t) = —j(t)v is a time-dependent drift force field, and a harmonic force is given by 
F(r) = -V'(r). The term F(t) is a fluctuating force with < F(t) >= and < r(ti)r(* 2 ) >= 
2DS(t 1 — t 2 ). From Eq. (5) we can find the dispersive behavior of a test particle, and that 
mathematical techique leads us to a more general result. 



(4) 



mr = f(v,t) - V'(r) + r(t) 



(5) 



In the three-dimensional case of an external harmonic potential V(r) 
solution of Eq. (5) is 



\ ar 2 , the formal 




(6) 



where Y(t) = exp[- /„ i(t')dt'] and j(t') = -y(t') + a. Then, one can show that a direct 
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calculation of the variance gives 



ft p—Zat 

o*(t) = De-**Y*(t)J Q y^dt (7) 

From Refs. 32 and 33, the quantity u 2 (t) is easily compared with that of other results. 
If 7(t) in terms of the time-dependent drift force field is given by a general form, we can 
discuss extentively the anomalous diffusive behavior of the tracer. Therefore, the Fokker- 
Planck equation associated with Eq. (5) is 

^P(f, t) = - V -[l(t)r + F(f)}P(r, t) + D V 2 P(r, t). (8) 

Especially, we will focus on two fractional Fokker-Planck equations: 

^P(r,t) = - V -[F(f,t)^P(f,t)] 

+ D V 2 ^P(r,t) (9) 



and 



^P(r,t) = - V .[F(r,t)P(r,t)]+ D V »P(r,t) (10) 



with | < r < 1 and 1 < v < 2, where r and v are, respectively, the time and the space 
scaling exponents. In Eq. (10), the Fourier-Laplace transform of the probability density 
function is described as 

P(k,u) = LF[P(r,t)Q(t)}, (11) 

where Q(t) is a step function, and F and L are, respectively, the Fourier and the Laplace 
transform operators. In the Fourier-Laplace transform domain, Eqs. (9) and (10) become 

[-iuj + D^i-iujf-^Pik.uj) = 
ik(—iu)) l ~ T 



(27T) 



d+l 



J dk'du'F(k- k',u-u')P(k',u') 

(12) 



and 

[-IU + D\k\ u ]P(k,u) = 



ik 



(27T) 



d+1 



y dk' du'F{k - k' ,u - J)P{k' ,u'). 

(13) 



Since A(fc,o;) and S(fc, o>) in the two cases where Eqs. (2) and (3) are, respectively, defined 

by 

A(jfe, cu) = — E(k, cu) = — (14) 

(—iuj) T + D\k\ 2 (-icu) T + D\k\ 2 

and 

A(fc, cu) = — , E(k, u) = (15) 

-iu + D\k\" -iu + D\k\» 

Eqs. (12) and (13) can be put into a unitary form: 



P(k,u) = A(k,u)- 



x J dk du G(k,k — k ,lo — to )P(k ,u ), 



(16) 



where 

G(k, k — k , cu — cu ) = S(fc, uj)k • F{k — k ,u — lu ). (17) 

For the special case of a tracer in the presence of a time-dependent force, the force-force 
correlation can be described as 

< F(f,t)F(r,t) >= b 2 (tt'fs(f-r)S(t-t) (18) 

and 

< F(k,u)F(k' ,u) >= 27rb 2 T(2{3 + l) 
x (-iu-iou'y m+1) 5(k 1 + k 2 ). (19) 

Hence, from the inverse Fourier-Laplace transform, the n-th moment of the tracer can be 
given by 
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where L~ x is the inverse Laplace-transform operator. In order to investigate the relaxation 
dynamics of the probability density toward the stationary solution, we now restrict our 
attention to the special class with the space scaling value £ = 1 and the time scaling value 
r = 2. We find the statistical behavior of a tracer from the moments of the probability 
density. From our fractional Fokker-Planck equation via Eqs. (12) — (19), the mean second 
moment for large t is 

From Eq. (21), < x 2 (t) > oc t with £ = \ for normal diffusion and < x 2 (t) > oc t v with r) > 1 
and £ > | for super diffusion. 

In conclusion, from a fractional Fokker-Planck equation, the tracer dipersion for an 
anomalous transport process has been discussed in the presence of a temporal power-law 
dependence of the drift and an external harmonic potential. The mean moments of the 
tracer were analyzed for the subdiffusive or superdiffusive behavior for several values of the 
exponent £ in the temporal power-law dependence of the drift. We think that for arbitrary 
time and space scaling exponents, it may not be easy to discuss the diffusive behaviors of the 
fractional Fokker-Planck equation. We also will present the results for non-Fickian diffusion 
for arbitrary time and space scaling exponents in other journals. We expect the detailed 
description of the anomalous behavior to be used to study extensions of numerical stratified 
models and fractal lattice models. 
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